Longitudinal evaluation of advanced glaucoma: ten year follow-up cohort study

This study focused on patients with advanced open-angle glaucoma (OAG) and aimed to identify key factors for monitoring them. We included 127 such patients who were followed for seven years or more, undergoing annual ophthalmic examinations. Glaucoma progression was defined as a deterioration in either structure or function. The progression rates and risk factors were evaluated. The patients were divided into upper- and lower-half subgroups based on the reduction in intraocular pressure (IOP) from the baseline. Over an 11-year period, glaucoma progression was detected in 59 eyes (46.5%). The rate of change in mean deviation (MD) was − 0.43 dB/year for the entire population; − 0.67 dB/year for progressors; and − 0.20 dB/year for non-progressors. Hypertension and disc hemorrhage (DH) were more common in progressors compared to non-progressors (45.8 vs. 23.5%, 11.9 vs. 1.5%; P = 0.008 and P = 0.016). Multivariate Cox’s proportional hazard model revealed that the presence of DH and a better baseline MD were associated with glaucoma progression. Additionally, patients with a higher percentage reduction in IOP (> 20.94%) had a lower risk of progression compared to those with less reduction. Inadequate IOP reduction, better baseline MD, presence of DH, and lower central corneal thickness were identified as risk factors for progression in advanced OAG patients.


Comparison of clinical characteristics of progressors and non-progressors
Of the 127 subjects, 59 eyes (46.5%) showed glaucomatous deterioration and were defined as progressors.Table 1 shows a comparison of the progressors' and non-progressors' clinical characteristics.In the progressors, the mean baseline RNFL and ganglion cell-inner plexiform layer (GCIPL) thicknesses were significantly better than those in the non-progressors (65.92 ± 8.65 vs. 61.34± 10.15 μm and 63.72 ± 8.04 vs. 60.83 ± 7.63 μm; P = 0.007 and P = 0.044, respectively).The mean deviation (MD) at baseline also differed between the two groups, with progressors presenting the better values (− 12.77 ± 4.17 vs. − 16.80 ± 5.56 dB; P < 0.001).The proportions of patients with underlying hypertension and DH also were higher in progressors (45.8 vs. 23.5% and 12.1 vs. 0%; P = 0.008 and P = 0.016, respectively).The results of the analysis of the rates of change of the structural and functional parameters are shown in Table 2.The rate of MD change was significantly faster for progressors than for nonprogressors (-0.67 ± 0.49 vs. -0.20 ± 0.25 dB/year; P < 0.001).

Risk factors associated with glaucoma progression
The results of the univariate and multivariate Cox's proportional hazard models for the effect of each variable on disease progression are detailed in Table 3.The univariate analysis showed that progression was significantly associated with presence of disc hemorrhage (DH) (P = 0.004) and baseline MD (P < 0.001).The stepwise multivariate Cox analysis included all of the parameters for which the P-value of the association with VF progression was ≤ 0.1 in the univariate analysis.Glaucoma progression was significantly associated with better baseline VF status (hazard ratios (HR) = 1.109; 95% confidence intervals (CI), 1.044-1.177;P = 0.001) and the presence of DH (HR = 2.352; 95% CI, 1.041-5.316;P = 0.040).The Kaplan-Meier survival analysis showed that, regardless of whether the progression criteria were structural or functional, only functional, or both, subjects with worse baseline MD had a higher cumulative non-progression probability than did those in the better baseline MD group (log rank test, P < 0.001, P < 0.001, and P = 0.012 for each progression criterion, respectively; Fig. 1 (A,B,D)).However, the analysis was inconsistent only for the structural progression criteria (log rank test, P = 0.110; Fig. 1C).In the group with better baseline MD, the lowest baseline MD was −13.55 dB.

Comparison of IOP reduction and glaucoma progression
Patients in the upper-half of the percentage-IOP reduction scale had a higher cumulative probability of nonprogression to glaucoma than those in the lower-half group, according to a Kaplan-Meier survival analysis (log rank test, P = 0.019; Fig. 2A).Supplementary Figure S1 displays Kaplan-Meier survival plots based on both progression criteria and only structural or functional progression criteria (all P's > 0.05).The upper-half group's lowest IOP reduction percentage was 20.94%.The upper-quartile percentage IOP reduction group (lowest IOP reduction of 31.0%) also revealed a higher cumulative probability of non-progression to glaucoma than did those in the lower reduction group, according to a Kaplan-Meier survival analysis (log rank test, P = 0.048; Fig. 2B).The Kaplan-Meier survival plots for the groups representative of structural or functional progression alone and both structural and functional progression also were examined.However, the percentage reduction of IOP did not show any significant differences for those progression criteria (all P's > 0.05, Supplementary Figure S2.A www.nature.com/scientificreports/multiple regression model revealed significant interactions of percentage-IOP reduction, central corneal thickness (CCT), and baseline standard automated perimetry (SAP) MD with the MD slope (β = 0.014, P = 0.003; β = 0.003, P = 0.029 and β = -0.029,P = 0.002, Table 4).

Comparison of clinical characteristics of participants with and without low vision
Nineteen ( 19) eyes (15.0%) of patients exhibited low vision at the most recent follow-up.When comparing the clinical characteristics between the patients presenting low vision at the final visit and those not, the baseline BCVA among patients with low vision was 0.24 ± 0.30 logMAR, which was worse than that among patients without low vision (0.04 ± 0.13 logMAR; P < 0.001), while baseline average GCIPL thickness was significantly lower (58.11± 7.94 vs. 62.94 ± 7.73 μm; P < 0.022, Supplementary Table S2).As for IOP fluctuation, it was significantly greater in the low-vision group (2.49 ± 1.12 vs. 2.03 ± 0.80; P = 0.034).

Risk factors associated with low vision
Medical history of hypertension, IOP fluctuation and baseline BCVA were found to be associated with visual impairment by the univariate Cox's proportional hazard model (Supplementary Table S3).In the multivariate Cox's proportional regression analysis employing a forward conditional method, poor baseline BCVA was significantly associated with presentation of low vision (HR = 8.915; 95% CI, 1.378-57.666;P = 0.022).

Representative cases
Representative cases of patients with advanced glaucoma are shown in Fig. 3 and 4. Figure 3 shows the longitudinal data for a male patient (age: 69) who demonstrated VF progression at the MD change of − 1.60 dB/year and the average RNFL thickness slope of −0.40 μm/year after 10 years of follow-up.Figure 4 displays the results for a male patient (age: 66) who demonstrated no signs of glaucoma progression during the 11-year follow-up period.

Discussion
This an 8-year period 9 .Kim et al. 10 retrospectively reviewed 87 eyes of medically treated advanced-normal tension glaucoma (NTG) patients with a baseline MD − 16.6 dB, and reported that 62.2% of them showed progression of VF damage over a period of 5.3 years.In their study, younger age and lower initial BCVA were found to be significant risk factors for loss of visual function.Recently, the Treatment of Advanced Glaucoma Study identified newly diagnosed advanced-glaucoma patients over the course of 36 months of follow-up at 27 centers 11 .
The trial found no evidence of any difference in vision-related quality of life between trabeculectomy or glaucoma drops 12 ; however, trabeculectomy significantly reduced the proportion of progressing eyes at 24 months 13 .Although clinical course is widely studied, to our knowledge, longitudinal investigations of 10 years' duration or longer are limited.In this study, we analyzed advanced glaucoma over the course of a mean follow-up period of 11 years.We believe that the long-term results of our study could offer deeper insight into end-stage glaucoma.The necessity of reducing IOP for treatment of glaucoma has been convincingly proven by previous investigations [14][15][16] .The AGIS results showed that patients who maintained an IOP of less than 18 mmHg over the course of six years had no VF progression 17 .According to our own findings (Fig. 2), the group that met the treatment goal of 20% IOP reduction was observed to have lesser glaucoma progression.However, the majority of the differences in Kaplan-Meier survival plots seem to have happened early, and thereafter the slopes parallelized.Long-term improvements in IOP control may lessen the impact of these early significant disparities.
In the present study, occurrence of DH during the study period significantly increased the risk of glaucoma progression.Usually, DH is thought to be frequently observed for the early to moderate stage of glaucoma and to show decreased incidence in optic discs with advanced cupping, especially in regions where the neuroretinal rim is absent 18,19 .The results of this study emphasize that the presence of DH is of vital importance, even for detection of progression of advanced glaucoma.Our results are consistent with previous studies in this regard [19][20][21][22][23] .In The worse baseline MD group (solid line) demonstrated a significantly greater cumulative probability of nonprogression than did the better baseline MD group (dotted line) for (A) the structural or functional progression criteria (P < 0.001), (B) the functionalprogression criteria (P < 0.001), or (D) both structural and functional progression criteria (P = 0.012).However, for the structural progression criteria, the cumulative probabilities of non-progression between the worse and better baseline MD groups were not significantly different (P = 0.110) (C).fact, development of DH may be linked to progression at all stages of glaucoma.Since progressors in the current study included a high percentage of patients with hypertension, some DH development may be associated with vascular factors.Furthermore, the inconsistency of the factors identified by the multiple regression model and Cox's proportional hazard models should be taken into account when considering the findings of our study.Different endpoint criteria might have contributed to the inconsistent results.
This study demonstrated that subjects with less severe VF defect at baseline were more likely to show glaucoma progression during follow-up.However, the relationship between baseline VF status and the possibility of glaucoma progression remains controversial.Several previous studies' results are consistent with ours in showing that better baseline VF status increases the risk of progression 9,24 .On the other hand, some other studies determined, contradictorily, that severe VF loss at baseline is a risk factor for progressive severe visual deterioration 7,[25][26][27] .A possible explanation for this discrepancy may be related to differences in study populations.All of our participants were classified as advanced glaucoma according to the Hodapp-Parrish-Anderson criteria, based on a mean baseline SAP MD of − 14.93 dB, which, compared with participants in other previous studies, did not represent severe damage 4 .Additionally, the "floor effect" of OCT measurements might be another reason contributing to the under-detection of disease progression for the individuals with severely damaged glaucoma.Further study will be necessary to elucidate the complex relationships between baseline severity of glaucoma and progression of VF damage.
Detecting progression of advanced glaucoma, which is an essential part of glaucoma management, is challenging due to the poor accuracy of structural thickness measurements 28 .The "floor effect" is defined as the point at which no further structural loss can be detected, and this, in fact, is thought to be a serious problem for monitoring of structural deterioration 28,29 .In this study, OCT was ineffective in detecting progression for those who had worse baseline MD (Fig. 1C).There have been attempts made to estimate SD-OCT-measured average floors for RNFL thickness.Bowd et al. 30 estimated the measurement floor for structural analysis based on variability in stable glaucomatous eyes using an arbitrary variability cutoff, and demonstrated that the average floors were 38 μm for circumpapillary RNFL thickness and 38 μm for GCIPL thickness.According to other previous reports, the estimated average RNFL thickness measurement floor ranged from 44.9 to 53.7 μm [31][32][33][34] .In the present study, the baseline average RNFL and GCIPL thicknesses averaged 63.48 and 62.20 μm, and the average RNFL and GCIPL thicknesses at final visit were measured to 57.99 and 57.79 μm, respectively (Supplementary Table S1).The mean RNFL and GCIPL thicknesses tended to be greater than the estimated measurements of previous studies, suggesting that the eyes enrolled for the present study had not reached the floor.This might have been due to the fact that these glaucomatous eyes were less advanced at baseline than were those in the relevant previous investigations.
This study has several limitations which should be considered when interpreting our results.First, the medical records reviewed were from a single tertiary medical center, and so selection bias might have been operative.The inclusion of 11 patients (8.7%) who had undergone glaucoma surgery could have biased our results.Our results may also have been affected by the exclusion of 8 patients (6.3%) who did not have 5 reliable results from VF tests.Second, all of the participants in this study were Koreans.Racial differences may have contributed to the study subjects' comparatively younger age of glaucoma onset.Longitudinal studies including non-East Asian patients can afford further clarification and more comprehensive knowledge of advanced glaucoma.Lastly, all of the study participants were categorized as having advanced glaucoma, with a mean baseline SAP MD of − 14.13 dB.Consequently, there appears to be insufficient data available on patients with even more severe advanced glaucoma.
In conclusion, this study presents data on eyes with advanced glaucoma that had been followed during a mean period of 11 years.Our findings suggest that development of DH, better baseline VF status and inadequate IOP control are significant risk factors for advanced-OAG progression.This should be considered in monitoring patients with advanced glaucoma.

Methods
This study included subjects from an ongoing study's OAG and NTG cohort at the Glaucoma Clinic of Seoul National University Hospital, and the study protocol was approved and informed consents were waived due to the retrospective nature of this study by the Institutional Review Board (IRB) of Seoul National University Hospital (IRB no.2304-077-1422).All of the investigations adhered to the tenets of the Declaration of Helsinki.www.nature.com/scientificreports/

Study participants
In this longitudinal cohort study, we reviewed the medical records of patients who had been diagnosed with advanced OAG and followed up for more than 7 years.Advanced glaucoma is classified according to the 'severe' category of VF loss based on the Hodapp-Parrish-Anderson criteria 35 : (1) MD worse than − 12 dB, (2) at least 50% of the points depressed at P-value < 5 or 20% of the points at P-value < 1% in the pattern deviation plot, (3) at least 1 point in the central 5 degrees with a sensitivity of 0 dB, or (4) points within the central 5 degrees with sensitivity < 15 dB in both hemifields.We initially selected 166 patients who had visited Seoul National University Glaucoma Clinic from August 2001 to September 2015 and been diagnosed with advanced OAG.Of those, patients who met the following criteria were enrolled: (1) treatment and a follow-up period of at least 7 years after the diagnosis of advanced OAG and (2) at least 5 reliable SAP results acquired with the Humphrey Field Analyzer (Carl Zeiss Meditec, Dublin, CA, USA) with the 24-2 test pattern using the Swedish Interactive Thresholding Algorithm (SITA) standard test protocol.VF results with fixation losses < 20%, false negatives ≤ 15%, and false positives ≤ 15% were included as reliable VFs.The exclusion criteria were as follows: (1) any history of ophthalmic surgery (except uncomplicated cataract surgery and glaucoma surgery) or retinal laser photocoagulation; (2) occurrence of ophthalmic disease other than glaucoma that might affect visual function; (3) presence of pseudoexfoliation or pigment dispersion syndrome; (4) possibility of secondary glaucoma related to uveitis, and (5) poor image quality.A total of 127 subjects were included in the evaluation.If both eyes were eligible for the investigation, one eye was selected by computerized random selection.
Based on the most acknowledged and tolerated medical therapy, it was the goal for every patient to reduce IOP by at least 20% and to less than 18 mmHg.In cases where this was not accomplished, surgical interventions were combined with medication therapy to lower IOP.

Clinical and ophthalmic evaluations
All of the study participants underwent complete ophthalmic examinations including measurement of BCVA, refractive error, slit-lamp biomicroscopy, IOP measurements using the Goldmann applanation tonometer (Haag-Streit, Konig, Switzerland), gonioscopy, dilated funduscopic examination, CCT measurement (Pocket II Pachymeter Echo Graph; Quantel Medical, Clermont-Ferrand, France), axial length (AXL) measurement (Axis II PR; Quantel Medical, Inc., Bozeman, MT, USA), stereoscopic optic disc photography, red-free RNFL photography (TRC-50IX; Topcon Corporation, Tokyo, Japan), as well as optic nerve head and macular imaging by Cirrus SD-OCT (Carl Zeiss Meditec), and SAP 24-2 testing (Humphrey Field Analyzer; Carl Zeiss Meditec) according to the SITA standard strategy.SD-OCT and SAP 24-2 testing were performed annually.Fluctuation of IOP was defined as a standard deviation of IOP during the follow-up period, and percentage reduction of IOP was determined using the following formula 100 × (baseline IOP-mean IOP during follow-up)/baseline IOP.The World Health Organization (WHO) classification was used to determine the level of visual impairment as measured by visual acuity.According to the WHO standards, blindness was classified as BCVA worse than 3/60 or constriction of the central VF to less than 10 degrees and severe visual impairment as BCVA between 6/60 and 3/60.In this study, participants who had BCVA worse than 6/60 were defined as low vision based on VA, and were designated as being in the low-vision group.

Assessment of optic disc hemorrhage
DH was defined as an isolated hemorrhage located in the disc tissue, or in the peripapillary retina connected to the disc rim, that is not associated with optic disc edema, papillitis, diabetic retinopathy, central or branch retinal vein occlusion, or any other retinal disease 36 .Throughout the study, occurrence of DH was evaluated at every visit by two glaucoma specialists (Y.I.S. and J.W.J.) based on stereoscopic optic disc photography and red-free RNFL photography.Discrepancies between the two observers' findings were resolved by discussion and consensus.

Figure 1 .
Figure 1.Kaplan-Meier survival plot of glaucoma progression of study population as stratified by baseline mean deviation (MD).The better baseline MD group included those with baseline MD greater than -13.55 dB.The worse baseline MD group (solid line) demonstrated a significantly greater cumulative probability of nonprogression than did the better baseline MD group (dotted line) for (A) the structural or functional progression criteria (P < 0.001), (B) the functionalprogression criteria (P < 0.001), or (D) both structural and functional progression criteria (P = 0.012).However, for the structural progression criteria, the cumulative probabilities of non-progression between the worse and better baseline MD groups were not significantly different (P = 0.110) (C).

Figure 2 .
Figure 2. Kaplan-Meier survival plot of glaucoma progression of study population as stratified by percentage reduction of intraocular pressure (IOP).The upper-half group included those with percentage reduction of IOP > 20.9%, and the upper-quartile group included those with percentage reduction of IOP > 31.0%.The cumulative probabilities of non-progression in the upper-half group (solid line) and the lower-half group (dotted line) were significantly different (P = 0.019) (A), and the cumulative probabilities of non-progression in the upper-quartile group (solid line) and the lower-reduction group (dotted line) were significantly different as well (P = 0.048) (B).

Figure 3 .
Figure 3. Representative case of progression of patient with advanced glaucoma.After a follow-up of 10 years, a male patient (age: 69) showed VF progression of advanced glaucoma with an MD change of − 1.60 dB/year and an average RNFL thickness change of -0.40 μm/year.

Figure 4 .
Figure 4. Representative case of non-progression of patient with advanced glaucoma.After a follow-up 11 years, a male patient (age: 66) manifested no signs of glaucoma progression.The mean MD change was -0.30 dB/year and the average RNFL thickness change was -0.16 μm/year.

Table 2 .
Rates of Change of Structural and Functional Parameters.Values are mean ± standard deviation.

Table 3 .
Cox's proportional hazard models for progression of advanced glaucoma.Factors with P < 0.1 in the univariate analysis were included in the multivariate analysis.IOP = intraocular pressure; SAP = standard automated perimetry; MD = mean deviation; VFI = visual field index; OCT = optical coherence tomography; RNFL = retinal nerve fiber layer; GCIPL = ganglion cell-inner plexiform layer.
longitudinal cohort study investigated the clinical course of advanced-open-angle glaucoma (OAG) patients and assessed the progression rate and risk factors of glaucoma progression for them.During a follow-up period of 10 years or more, 46.5% of patients showed glaucoma progression with a mean MD change of − 0.67 ± 0.49 dB/ year.Better baseline VF status, presence of DH, inadequate IOP lowering, and lower CCT were factors associated with glaucoma progression.There have been previous attempts to demonstrate the clinical course and risk factors for advanced glaucoma.The Advanced Glaucoma Intervention Study (AGIS), the first multicenter randomized clinical trial concerning advanced glaucoma 8 , demonstrated that IOP-lowering extent was associated with progression of VF damage over Vol.:(0123456789) Scientific Reports | (2024) 14:476 | https://doi.org/10.1038/s41598-023-50512-7

Table 4 .
Multiple regression analysis for slope of mean deviation.Factors with P < 0.1 in the univariate analysis were included in the multivariate analysis.IOP = intraocular pressure; SAP = standard automated perimetry; MD = mean deviation; VFI = visual field index; OCT = optical coherence tomography; RNFL = retinal nerve fiber layer; GCIPL = ganglion cell-inner plexiform layer.